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and 



By the first-principles electronic structure calculations, we find that the ground state of ternary 
iron selenides AFe2Se2 (A=K, Cs, or Tl) is in a bi-collinear antiferromagnetic order, in which the 
Fe local moments 2.8/jb) align ferromagnetically along a diagonal direction and antiferromag- 
netically along the other diagonal direction on the Fe-Fe square lattice. This bi-collinear antiferro- 
magnetic order results from the interplay among the nearest, the next nearest, and the next next 
nearest neighbor superexchange interactions, mediated by Se 4p-orbitals. 

PACS numbers: 74.25.Jb, 71.18.+y, 74.70.-b, 74.25.Ha, 71.20.-b 



The discovery of high transition temperature Tc su- 
perconductivity in LaFeAsO by partial substitution of O 
with F atoms [1] stimulates the intense studies on the iron 
pnictides. There have been four types of iron-based com- 
pounds reported to show superconductivity after dop- 
ing or under high pressures, i.e. 1111-type i^eFeAsO 
{Re = rare earth) [ij, 122-type 5Fe2As2 (5=Ba, Sr, or 
Ca) [2|, Ill-type AFeAs {A = alkah metal) [3], and 11- 
type tetragonal a-FeSe(Te) [1]. It was very recently re- 
ported that the superconductivity was observed at about 
30 K in a new FeSe-layer compound Ko.8Fe2Se2 [sf, 
formed by intercalating potassium (K) atoms between 
FeSe layers. Soon after, the superconductivity was also 
found in the Cs- intercalated compound Cso.8(FeSeo.98)2 
[gI] and Tl-intercalated compound TlFea^Se2[7]. Although 
these the compounds take the ThCr2Si2 type struc- 
ture with PA/nmm symmetry isostructural with 122- 
type 5Fe2As2, they can be considered as a new type of 
iron-based superconductors since they are chalcogenides 
rather than pnictides. 

For the parent compounds AFe2Se2 (A=K, Cs, or 
Tl) of these new superconductors, the intercalated alkali 
metal or Tl atoms will directly dope much more electrons 
into the FeSe layers than for the iron pnictides. It is thus 
expected that the electronic and magnetic structures of 
these new compounds are very likely different from the 
ones of those iron pnictides. On the other hand, in or- 
der to investigate the mechanism of superconductivity in 
these materials, one needs to first understand the elec- 
tronic and magnetic structures of the parent compounds 
AFe2Se2. 

In this paper, we report the theoretical result on 
the electronic and magnetic structures of iron selenides 
AFe2Se2 (A=K, Cs, or Tl) obtained from the first- 
principles electronic structure calculations. We find that 
the compounds AFe2Se2 are antiferromagnetic semimet- 
als with a bi-collinear antiferromagnetic order in the 
ground states, resulting from the strong nearest (Ji), 
next-nearest (J2), and next-next-nearest (J3) neighbor 
superexchange interactions in these materials. A small 
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FIG. 1: (Color online) CsFe2Se2 with the ZrCuSiAs-type 
structure: (a) a tetragonal unit cell containing two formula 
units; (b) schematic top view of the FeSe layer. The large 
dashed square is a a x 2a unit cell. The Fe spins in the bi- 
collinear antiferromagnetic order are shown by the red arrows. 



monoclinic lattice distortion due to spin-lattice coupling 
was further found in these compounds, similar to a-FeTe 
[sl, m. Here the bi-collinear antiferromagnetic (AFM) 
order means that the Fe moments align ferromagneti- 
cally along a diagonal direction and antiferromagnetically 
along the other diagonal direction on the Fe-Fe square 
lattice. In other words, if the Fe-Fe square lattice is di- 
vided into two square sublattices, the Fe moments on 
each sublattice take their own collinear AFM order. The 
bi-collinear antiferromagnetic order was first predicted 
for a-FeTe by our previous work [H and confirmed by 
the later neutron scattering experiment [9]. 

Although AFe2Se2 can be considered as a tetragonal 
crystal with two formula units included in the corre- 
sponding unit cell as shown in Figdja), its primitive unit 
cell is constructed by considering AFe2Se2 as a triclinic 
crystal, in which only one formula unit cell is included. In 
the nonmagnetic calculations, we adopted the primitive 
cell as a calculation cell. In the calculation of the elec- 
tronic and magnetic structures of the ferromagnetic and 
the square antiferromagnetic Neel states, the a x a FeSe 
cell is taken as the base cell. In the calculation of the 
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FIG. 2: (Color online) CsFe2Se2 in the nonmagnetic state: 
(a) energy band structure; (b) the Fermi surface sheets due 
to the band crossing the Fermi energy marked by the blue 
line in the left figure; (c) the Fermi surface sheet due to the 
band crossing the Fermi energy marked by the red line in the 
left figure. 



collinear and bi-collinear antiferromagnetic states, the 
unit cells are doubled and the base cells are the >/2a x \f2a 
and a X 2a (Fig. [11(b)) FeSe cells, respectively. 

In our calculations the plane wave basis method was 
used [lo|- We adopted the generalized gradient approx- 
imation (GGA) of Perdew-Burke-Ernzerhof pjj for the 
exchange-correlation potentials. The ultrasoft pseudopo- 
tentials [lH were used to model the electron-ion interac- 
tions. After the full convergence test, the kinetic energy 
cut-off and the charge density cut-off of the plane wave 
basis were chosen to be 800 eV and 6400 eV, respectively. 
The Gaussian broadening technique was used and a mesh 
ofl8xl8x9 k-points were sampled for the Brillouin- 
zone integration. In the calculations, the experimental 
tetragonal lattice parameters were adopted P, IgI] , and the 
internal atomic coordinates within the cell were deter- 
mined by the energy minimization. Actually, the lattice 
parameters optimized by the energy minimization in our 
magnetic calculations are found in excellent agreement 
with the experimental ones. 

In Fig. [2] we plot the electronic band structure and 
Fermi surface of CsFe2Se2 in the nonmagnetic state, simi- 
lar to the ones just reported for KFe2Se2 ^3|. As we see, 
there are three Fermi surface sheets, among which the 
two cylinder-like ones are located around the corners of 
the Brillouin zone, and another one is of small pocket. 
The Fermi surface is all of electron- type. The volume 
enclosed by these Fermi sheets give 1.008 electrons per 
formula cell. i.e. about 0.84 x 10^^/cm^ . The den- 
sity of states (DOS) at the Fermi energy is about 2.27 
states per eV per formula unit. The corresponding elec- 
tronic specific heat coefficient and Pauli susceptibility are 
7 = 5.35 mJ/ (K'^^mol) and x = 0-92 x 10~^ m^/mol, re- 
spectively. Our calculations for the nonmagnetic states 
exclude any possible structural distortions. This sug- 
gests that any possible structural distortion happening 
in CsFe2Se2 would be driven through spin-phonon inter- 
actions. For KFe2Se2, the almost same results are found 
by our calculations. 



In order to explore the magnetic structure of CsFe2Se2, 
we have calculated four different possible magnetic states 
with ferromagnetic, square Neel AFM, collinear AFM, 
and bi-collinear AFM orders, respectively. If the en- 
ergy of the nonmagnetic state is set to zero, we find 
that the energies of the ferromagnetic {Efe)-, square Neel 
AFM [Eafm), collinear AFM [Ecol), and bi-collinear 
AFM states {Ebi) are (-0.2122, -0.1651, -0.2677, -0.2752) 
eV/Fe for CsFe2Se2 and (-0.2437, -0.1544, -0.2613, - 
0.2884) eV/Fe for KFe2Se2, respectively. Thus the 
ground states of both CsFe2Se2 and KFe2Se2 are in the 
bi-collinear antiferromagnetic order, similar to the one 
of Qf-FeTe [11 . The magnetic moment around each Fe 
atom is found to be about 2.4 ~ 3.0 /i^, varying weakly 
in the above four magnet icallyordered states, similar as 
in LaFeAsO and BaFe2As2 3 [l5[. Since these local 
moments are embedded in the environment of itinerant 
electrons, the moment of Fe ions is thus fluctuating. 

To quantify the magnetic interactions, we assume that 
the energy differences between these magnetic orderings 
are predominantly contributed from the exchange inter- 
actions between the Fe moments with spin 5, which can 
be effectively modeled by the following frustrated Heisen- 
berg model with the nearest, next-nearest, and next next 
nearest neighbor couplings Ji, J2, and J3, 



{{{ij))) 



whereas (zj), <C ij ^ and {{{ij))) denote the sum- 
mation over the nearest, next-nearest, and next-next- 
nearest neighbors, respectively. This model may miss 
certain contributions from itinerant electrons, however, 
we believe that it captures the substantial physics on the 
magnetic structures. From the above calculated energy 
data, we find that for CsFe2Se2 Ji = -11.78 meV/S'^, 
J2 = 19.75 meV/S^, and J3 = 11.75 meV/S'^, while for 
KFe2Se2 Ji = -22.3 meV/S'^, J2 = 15.6 meV/S'^, and 
J3 = 14.6 meV/S'^ (The detailed calculation is referred 
to Appendix of Ref. [l3 ). Notice that the ferromagnetic 
states in both the compounds are lower in energy than 
the square Neel AFM state, which results in the negative 
Ji here. 

It is known that when J3 > J2/2 and J2 > Ji/2, 
the bi-collinear AFM state is lower in energy than the 
collinear AFM state for a frustrated J1-J2-J3 Heisenberg 
model. This is in well agreement with the derived Ji, 
J2, and J3 for both CsFe2Se2 and KFe2Se2. Accord- 
ingly, by J1-J2-J3 Heisenberg model we may understand 
the complex magnetic structures displayed in both the 
compounds. On the other hand, the classical study of 
J1-J2-J3 Heisenberg model ll,[l7| shows that there may 
be an incommensurate AFM spin order easily developed 
in such J1-J2-J3 system. However, it can be shown that a 
slight structural distortion making Ji different along dif- 
ferent direction will eliminate such an incommensurate 
spin order [ITJ. 
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FIG. 3: (Color online) (a) Electronic band structure of 
CsFe2Se2 in the bi-collinear-ordered antiferromagnetic state. 
The Fermi energy is set to zero, (b) Brillouin zone, (c) 
Fermi surface. Notice that r-N(r-N^) corresponds the 
parallel(antiparallel)-aligned moment line. 
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FIG. 4: (Color online) Calculated total and projected density 
of states at the five Fe-Sd orbitals around one of the four Fe 
atoms in the bi-collinear antiferromagnetic state unit cell of 
CsFe2Se2. The Fermi energy is set to zero. 



It is expected that there would be a further lattice 
distortion considering possible spin-phonon interactions. 
Similarly to spin-Peierls distortion, the lattice constant 
slightly expands along spin anti-parallel alignment to 
lower AFM energy and/or slightly contracts along spin- 
parallel alignment to lower further ferromagnetic energy. 
Indeed such small structural distortions are found for 
both CsFe2Se2 and KFe2Se2 with an extra energy gain 
of ~4 meV/Fe. As a result, the crystal unit cell of the 
compounds on FeSe layer deforms from a square to a 
rectangle (rhombus). Such small lattice distortions af- 
fect weakly the electronic band structures and the Fe 
moments. 

Our calculations also show that for both CsFe2Se2 and 
KFe2Se2 the Fe magnetic moments between the nearest 
neighbor layers FeSe prefer the parallel alignment with a 
small energy gain of about 2 meV/Fe in comparison with 
the anti-parallel alignment. This is different from the iron 
pnictides [H, [l5| . It is thus very likely that here the 
magnetic phase transition would happen simultaneously 
with the structural transition. Overall the magnetic or- 
der vector in these two compounds is thus (^, 0, 0). The 
corresponding magnetic Bragg peaks are (1,0,0). 

Fig. [3] shows the electronic structure of CsFe2Se2 in 
the bi-collinear AFM state, and the similar finding is 
for KFe2Se2. There are three bands crossing the Fermi 
level which form four sheets of the Fermi surface. The 
Fermi surface contains two hole-type sheets parallel to 
the plane F-Z-R'-N', and two electron- type cylinders 
parallel to the plane A-M-N-R. From the volumes en- 
closed by these Fermi surface sheets, we find that the 
electron (hole) carrier density is 0.052 electron/formula 
cell (0.061 hole/formula cell), namely, 0.88 x lO^V^^^ 
(1.00 X lO'^^/cm^). The density of states (DOS) at the 
Fermi level E^? is 6.24 states per eV per formula cell. The 
corresponding electronic specific heat coefficient 7 = 7.34 
mJ/{K'^ * mol). 

By projecting the density of states onto the five 3d 
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FIG. 5: (Color online) Total and orbital-resolved partial den- 
sity of states (spin- up part) per four formula cells of CsFe2Se2 
in the bi-collinear antiferromagnetic state. The Fermi energy 
is set to zero. 



orbitals of Fe of CsFe2Se2 in the bi-collinear AFM state 
(see Fig. [4]), we find that the five up-spin orbitals are 
almost filled and the five down-spin orbitals are nearly 
uniformly filled by half. This indicates that the crystal 
field splitting imposed by Se atoms is very small. As 
the Hund rule coupling is strong, this would lead to a 
large magnetic moment formed around each Fe atom, as 
found in our calculations. The formation of Fe magnetic 
moments is thus due to the Hund's rule coupling, which is 
a universal feature found in all the iron-pnictides 13, [HI . 

Further inspection of the charge distribution in real 
space shows that there is the covalence bond formed be- 
tween the two nearest neighbor Fe and Se atoms, which is 
responsible for the next nearest neighbor Fe-Fe superex- 
change coupling J2 . This superexchange is antiferromag- 
netic because the intermediated state associated with the 
hopping bridged by Se atom is a spin singlet. Thus we 
find the Se-bridged superexchange antiferromagnetic in- 
teraction between the next nearest neighbor Fe-Fe atoms 
in CsFe2Se2 and KFe2Se2, which is again similar to the 
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FIG. 6: (Color online) Electronic band structure of TlFe2Se2 
in the bi-collinear-ordered antiferromagnetic state. The Bril- 
louin zone is shown in Fig. [3jb). Notice that r-N(r-N') 
corresponds the parallel(antiparallel)-aligned moment line. 



ones in all the iron pnictides [IJ, ll5|- Nevertheless, there 
is a substantial difference between CsFe2Se2 (KFe2Se2) 
and the iron pnictides, regarding the next next nearest 
neighbor Fe-Fe exchange interaction J3 . 

Fig. [5] shows that the band formed by Se 4p orbitals 
is partially filled at the Fermi energy in CsFe2Se2 and 
KFe2Se2. So there are itinerant electrons at the Fermi 
energy involved in mediating the exchange interactions in 
CsFe2Se2 and KFe2Se2. This may explain why the next 
next nearest neighbor superexhange coupling J3 is large 
for CsFe2Se2 and KFe2Se2 due to an RKKY-like mecha- 
nism. In contrast, the band formed by As 4p orbitals is 
gapped at the Fermi energy in the iron pnictides 3, HI , 
correspondingly the next next nearest neighbor superex- 
change coupling J3 is nearly zero. 

A physical picture suggested in our study is thus that 
the Fe moments in CsFe2 and KFe2Se2 are also partially 
mediated by partially delocalized Se 4p-band and the ori- 
gin of J3 exchange coupling may be well induced through 
an RKKY-like mechanism besides the Se-bridged su- 
perexchange antiferromagnetic interaction J2. Mean- 
while, the exchange interactions J2 and J3 are in strong 
frustration competition, which is sensitively modulated 
by the Fe-Se bonding and related local environment, as 
we see from a-FeSe (J3 < /2/2) [8] to CsFe2Se2 and 
KFe2Se2 here (J3 > I2/2). By our calculation, it is clear 
that the magnetic order should be along the diagonal 
direction of the Fe-Fe square lattice for CsFe2Se2 and 
KFe2Se2 since J3 is larger than J2/2. However, it is also 
very likely that defects, vacancies, or other impurities 
would drive the bi-collinear antiferromagnetic order into 
an incommensurate spin order along the diagonal direc- 
tion in CsFe2Se2 and KFe2Se2, similar to the ones found 
in a-FeTe [18]. 

Although Tl is not an alkali metal element, our cal- 
culations show that the electronic and magnetic struc- 
tures of TlFe2Se2 are overall quite similar to the ones of 



CsFe2Se2 and KFe2Se2. The electronic band structure 
and the Fermi surface of TlFe2Se2 in the nonmagnetic 
state are found by our calculations to be the same as the 
ones reported for TlFe2Se2 in Ref. [l9|. However, we find 
that the ground state of TlFe2Se2 is in the bi-collinear 
AFM order with a magnetic moment of ~ 2.7 fiB around 
each Fe atom rather than the square Neel AFM order 
reported in Ref. 19. The corresponding electronic band 
structure is shown in Fig. [6l 

In conclusion, we have presented the results of the 
electronic band structure and magnetic properties of 
AFe2Se2 (A=K, Cs, or Tl) based on the first-principles 
electronic structure calculations, our studies show that 
the ground state of AFe2Se2 is a quasi-2-dimensional bi- 
collinear antiferromagnetic semimetal with a magnetic 
moment of ~ 2.8/ib around each Fe atom. This bi- 
collinear antiferromagnetic state can be understood by 
the Hesienberg model with J1-J2-J3 superexchange in- 
teractions. 
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